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TURBULENCE  IN  THE  WAKE  OF  A  THIN  AIRFOIL  AT  LOW  SPEEDS 
By  George  S.  Campbell 

SUMMARY 


Experiments  have  been  made  to  determine  the  nature  of  turbulence 
in  the  wake  of  a  two-dimensional  airfoil  at  low  speeds.  The  experi- 
ments were  motivated  by  the  need  for  data  which  can  be  used  for  analysis 
of  the  tail- buffeting  problem  in  aircraft  design.  Turbulent  intensity 
and  power  spectra  of  the  velocity  fluctuations  were  measiired  at  a 

Reynolds  number  of  1.6  x  10^  for  several  angles  of  attack.   Total-head 
measiirements  were  also  obtained  in  an  attempt  to  relate  steady  and 
fluctuating  wake  properties. 

Mean-square  downwash  was  found  to  have  nearly  the  same  dependence 
on  vertical  position  in  the  wake  as  that  shown  by  total-head  loss. 
For  this  particular  wing,  turbulent  intensity,  integrated  across  the 
wake,  increased  roughly  as  the  3/2  power  of  the  drag  coefficient. 

Power- spectrum  measurements  indicated  a  decrease  in  frequency  as 
wing  angle  of  attack  was  increased.  The  average  frequency  in  the  wake 
was  proportional  to  the  ratio  of  mean  wake  velocity  to  wake  width. 


INTRODUCTIC»J 


The  study  of  velocity  fluctuations  in  the  wake  of  an  airfoil  is 
of  interest  in  connection  with  the  buffeting  problem.  In  one  forni  of 
buffeting,  referred  to  as  tail  buffeting,  the  horizontal  tail  of  the 
aircraft  is  immersed  in  the  wake  of  the  wing  for  certain  angles  of 
attack.   If  the  wake  contains  sufficient  turbulent  energy,  as  it  may 
due  to  separated  flow  over  the  wing,  the  aircraft  may  experience  un- 
desirable load  fluctuations.  It  might  be  mentioned  that  buffeting  is 
not  always  undesirable,  as  it  is  often  useful  as  an  advance  warning 
against  dangerous  stability  deficiencies. 

A  recent  book  on  aeroelasticity  (ref.  1)  has  summarized  the  re- 
siilts  of  some  investigations  of  tail  buffeting.  From  this  summary, 
it  appears  that  there  has  been  considerable  work  done  on  buffeting 
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for  high  angles  of  attack,  but  that  there  is  little  information  avail- 
able on  the  possibility  of  biiffeting  at  lower  angles.  The  use  of 
increasingly  thinner  wing  sections  introduces  the  likelihood  of  leading- 
edge  separation  at  low  angles  of  attack  (ref.  2).  Although  this  type  of 
separation  usually  disappears  at  transonic  Mach  numbers,  it  is  generally 
replaced  by  separation  behind  the  main  compression  shock.   In  either 
case,  the  occurrence  of  such  separated  flows  at  low  angles  of  attack 
would  be  expected  to  cause  buffeting  difficulties. 

A  usefiil  approach  for  the  study  of  most  of  the  varioxis  types  of 
tail  buffeting  has  been  provided  in  reference  3.   In  this  paper,  buffet- 
ing was  considered  to  be  the  response  of  a  linear  system  to  a  random 
forcing  function.  In  order  to  predict  the  response  of  an  elastic  stinic- 
ture  to  turbulent  velocity  fluctuations,  one  must  know  the  power  spectr\im 
of  the  turbulence,  the  aerodynamic  impedance  of  the  lifting  surface,  and 
the  frequency  response  of  the  structure.  Nonstationary  wing  theory  may 
be  used  to  estimate  the  aerodynamic  impedance  functions  of  various  lift- 
ing surfaces  .  The  determination  of  the  frequency-response  character- 
istics of  aircraft  structures  has  received  sufficient  attention  over  a 
period  of  years  so  that  this  portion  of  the  buffeting  problem  may  be 
considered  solved.  The  remaining  quantity  necessary  for  the  analysis 
of  tail  buffeting  is  the  power  spectriim  of  turbulence  in  the  wake  of  the 
forward  lifting  surface.   Up  to  the  present  time,  it  has  been  necessary 
to  assume  that  the  turbulence  is  isotropic  or  that  it  possesses  a  dis- 
crete frequency  corresponding  to  a  Karman  vortex  street.  However,  it 
is  evident  that  measurements  of  the  power  spectra  of  turbulence  behind 
wings  are  necessary  before  the  general  statistical  approach  to  buffeting 
can  be  fully  utilized. 

The  primary  purpose  of  the  present  investigation  has  been  the 
measurement  at  low  speeds  of  the  power  spectrum  of  turbulence  behind  a 
thin,  two-dimensional  airfoil.  The  tests  were  made  by  means  of  the 
usual  hot-wire  techniques  for  measuring  velocity  fluctuations  in  turbu- 
lent flow  (ref.  5).  Angle  of  attack  was  varied  from  0  to  22  degrees, 

and  the  Reynolds  number  for  most  of  the  tests  was  about  1.6  x  10  .   It 
might  be  expected  that  the  results  could  be  applied  at  higher  Reynolds 
numbers  provided  the  wing  exhibits  the  type  of  leading-edge  separation 
shown  by  the  airfoil  in  the  present  tests.  Total-head  surveys  and 
surface-tuft  studies  were  made  in  order  to  determine  the  static  character- 
istics of  the  wing  and  its  wake.  Wherever  possible,  an  attempt  has  been 
made  to  relate  the  power  spectrum  of  turbulence  to  the  mean  properties 
of  the  flow. 


A  recent  experimental  investigation  (ref.  h)  confinns,  in  general, 
the  applicability  of  classical  nonstationary  wing  theory  when  the  veloc- 
ity fluctuations  are  of  a  randcm  nature. 


NACA  TM  ll|27 


This  research  was  conducted  at  the  California  Institute  of  Tech- 
nology under  the  guidance  of  H.  W.  Liepraann  while  the  author  was  a 
holder  of  a  Howard  Hughes  Fellowship. 


APPARATUS  AND  TESTS 


Wind  Tunnel 


The  experiments  of  this  investigation  were  made  in  the  GALCIT 
20-  by  20-inch  low  turbulence  wind  txmnel;  a  diagram  of  the  tunnel  is 
presented  in  reference  6.  The  turbulence  level  of  the  tunnel  is  about 
0.03  percent. 

Most  of  the  present  tests  were  conducted  at  velocities  of  about 
50  feet  per  second.   Free-stream  velocity  was,  in  general,  measured 
using  a  pitot-static  tube  and  a  micromanometer.  At  lower  speeds, 
however,  it  was  found  convenient  to  measiure  velocity  by  observing  on 
an  oscilloscope  the  frequency  of  vortices  shed  from  a  reference  cylin- 
der. A  description  of  this  procedure  is  given  in  reference  7. 

The  traverse  mechanism  for  positioning  hot-wire  and  total-head 
probes  is  the  same  as  the  one  described  in  reference  7. 


Airfoil 


A  two-dimensional  airfoil  of  6-inch  chord  spanned  the  test  section. 
The  airfoil  section  was  a  l/8-inch  thick  flat  plate  with  a  blunt  trail- 
ing edge  and  a  rounded  nose.  The  leading-edge  radius  was  approximately 
half  the  airfoil  thickness.  The  wing  was  made  of  steel  and  was  tested 
in  a  smooth  condition. 

All  measurements  were  made  behind  the  center  section  of  the  wing 
span.  Reynolds  number  for  most  of  the  tests  was  about  1.6  x  10  . 


Total-Head  Surveys 


A  total-head  survey  of  the  wake  was  made  one  chord  length  behind 
the  wing  trailing  edge  using  a  total-head  tube  of  0.108-inch  outer 
diameter.  Profile-drag  coefficient  was  obtained  by  integrating  the 
total-head  profile 
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This  formula  is  expected  to  give  quite  satisfactory  accuracy  with  the 
survey  plane  one  chord  from  the  trailing  edge  (ref.  8). 


Hot-Wire  Apparatus 

The  usual  methods  of  hot-wire  anemometry  (ref.  $)  were  used  to 
measure  fluctuating  velocities  in  the  airfoil  wake.  Hot  wires  of  0.1 
mil  platinum  or  platinum-rhodixim  were  used  for  all  measurements.  For 

measurements  of  the  longitudinal  velocity  fluctuation,  u  ,  a  l/8-inch 
long  wire  was  soldered  to  the  tips  of  the  probe.  Transverse  velocity 
fluctuations  were  obtained  with  a  probe  having  two  wires  of  0.2  inch 
length  inclined  at  angles  of  +  30  degrees  from  the  stream  direction. 

A  considerable  saving  in  testing  time  was  accomplished  by  omitting 
hot-wire  calibrations.  Reliable  relative  measurements  of  the  turbulent 
velocity  fluctuations  could  be  obtained  in  this  manner  provided  the 
hot-wire  apparatus  was  linear.   It  is  estimated  on  the  basis  of  data 
from  reference  9  that  the  downwash  probes  of  the  present  investigation 

should  be  linear  for  fluctuations  up  to  about  ±  2^  .  The  overall 
linearity  of  the  hot-wire  apparatus  was  checked  by  observation  of 
oscilloscope  traces  of  the  hot-wire  output.  The  presence  of  nonlin- 
earities  (such  as  those  which  could  be  introduced  by  excessively  high 
amplifier  gains)  generally  resulted  in  flat  tops  to  an  otherwise  random 
oscillograph  trace. 

The  frequency  response  of  the  amplifier  used  for  the  hot-wire 
measurements  was  flat  for  frequencies  of  l/2  to  25,000  cycles  per 
second.  The  post-amplifier  was  provided  with  resistance-capacitance 
compensation  which  was  adjusted  by  the  usual  square-wave  method.  In 
order  to  minimize  noise  a  10  kilocycle  cut-off  filter  was  used  in  the 
present  measurements. 


Electronic  Measurements 


Intensity  of  the  turbulent  velocity  fluctuations  was  determined 
by  means  of  an  average-reading  vacuum-tube  voltmeter.  Although  the 
voltmeter  was  calibrated  to  give  RMS  values  for  sinusoidal  voltage 
variations,  the  instrument  did  not  provide  true  RMS  intensities  for 
random  inputs.  On  the  basis  of  turbulence  measurements  with  a  statis- 
tical analyzer  (ref.  7),  the  difference  between  true  and  indicated 
RMS  is  expected  to  be  less  than  10  percent. 

Power-spectral  density  was  measured  with  a  Hewlett-Packard  model 
30OA  wave  analyzer.  A  constant  half -band  width  of  30  cycles  per  second, 
based  on  an  attenuation  of  UO  decibels,  was  used  for  all  measurements. 
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The  filtering  characteristics  of  the  wave  analyzer  are  shown  in  figure 
1  for  the  30  cycle  band  width.  The  frequency  range  of  the  analyzer  was 
about  30  -  16,000  cycles  per  second.  The  output  was  read  directly  from 
the  instrument's  voltmeter,  and  so  the  values  obtained  were  subject  to 
the  same  errors  for  random  inputs  as  were  the  intensity  measurements 
described  in  the  previous  paragraph. 


ACCURACY 


Fan  Frequency 


Of  the  factors  affecting  the  accuracy  of  measurements  of  power- 
spectral  density,  the  rotational  speed  of  the  wind-tunnel  fan  is  one  of 
the  most  important.  The  motor  and  fan  were  quite  well  isolated  from 
the  test  section  so  that  vibrations  were  not  transmitted  to  the  hot- 
wire mount.  However,  periodic  pressure  fluctuations  were  found  to  be 
transmitted  upstream,  and  so  a  test  was  made  to  determine  the  effect 
of  tunnel  fan  rotation  on  power-spectrum  measurements.  Downwash  veloc- 
ities are  of  primary  interest,  and  so  downwash  spectra  are  compared 
in  figure  2  for  two  fan  speeds.  The  same  free-stream  velocity  was 
maintained  by  introducing  a  screen  just  ahead  of  the  fan  for  one  of  the 
tests.  Without  the  screen,  the  peak  in  the  downwash  spectrum  occurs  at 
nearly  the  same  frequency  as  the  fan's  second  harmonic.  However,  the 
downwash  spectrum  has  practically  the  same  shape  when  the  fan  and  wing 
peaks  do  not  coincide,  i.e.,  when  the  screen  is  added.  Therefore,  it 
may  be  concluded  that  the  fan  is  not  the  cause  of  the  peak  in  the  down- 
wash  spectrxtm,  even  though  it  does  somewhat  modify  the  spectrum. 


Noise 


The  noise  introduced  by  the  measuring  apparatus  was  generally 
small  and  was  subtracted  from  the  total  output.  Near  the  center  of  the 
wake,  the  ratio  of  signal  to  noise  was  greater  than  100  to  1  for  RMS 
turbulent  intensity.  Measurements  at  60  cycles  per  second  were  usually 
high  by  about  3  percent  due  to  pickup  of  line  frequency.  The  greatest 
luicertainty  in  the  spectrum  measurements  occurred  at  the  lowest  fre- 
quency, 30  cycles  per  second,  due  to  the  characteristics  of  the  wave 
analyzer. 


Wall  Interference 


On  the  basis  of  an  analysis  by  Glauert  (ref.  10),  the  Strouhal 
number  of  the  vortex  street  obtained  at  an  angle  of  attack  of  21.7 
degrees  in  the  present  investigation  should  be  about  6  percent  higher 
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than  the  value  obtained  in  a  stream  of  unlimited  extent.  This  is  also 
the  magnitude  of  the  correction  necessary  to  bring  the  present  results 
into  agreement  with  those  of  Fage  and  Johansen  (ref.  11),  obtained  for 
the  same  size  model  in  a  much  larger  tunnel.  The  corrections  to 
Strouhal  number  at  lower  angles  of  attack  are  proportionately  smaller. 
It  appears  that  wall  corrections  for  power  spectra  and  turbulent  in- 
tensities are  not  presently  available.   In  order  to  have  consistency 
between  Strouhal  numbers  and  power  spectra,  all  data  are  therefore  pre- 
sented without  applying  any  corrections  for  wind-tunnel  wall  interference. 


Hot-Wire  Length  Corrections 


The  general  effect  of  the  finite  length  of  a  hot  wire  is  to  reduce 
the  voltage  across  the  wire  as  a  result  of  imperfect  correlation  of  the 
velocity  fluctuations  at  different  points  of  the  wire.  For  isotropic 
turbulence,  estimates  of  the  corrections  to  intensity  may  be  quite 
simply  obtained  (references  12  and  13,  for  example).   In  experimental 
measurements  of  the  spectrum  of  turbulence  behind  grids  (ref.  6),  cor- 
rections to  the  power-spectral  density  were  found  negligible  for  reason- 
able hot-wire  lengths.  Unfortunately,  the  estimation  of  length  correc- 
tions for  the  present  investigation  is  greatly  complicated  by  the  use 
of  a  two-wire  probe  in  a  field  which  bears  scant  resemblance  to  isotropic 
turbulence.   In  view  of  this  uncertainty,  length  corrections  have  not 
been  applied  to  the  data. 


RESULTS 


Wake  Energy  Distribution 


Two  measures  of  wake  energy  were  obtained  in  the  present  investiga- 
tion: total-head  loss  and  the  mean-square  of  the  downwash  fluctuations. 
As  total-head  measurements  are  considerably  simpler  than  hot-wire 
studies,  it  is  desirable  to  determine  whether  properties  of  the  turbu- 
lent flow  can  be  related  to  the  simpler  mean  flow. 

Presented  in  figures  3-7  are  the  wake  profiles  obtained  from 
total-head  and  hot-wire  surveys  at  a  distance  of  one  chord  length  from 
the  wing  trailing  edge.  The  profiles  have  been  normalized  to  unit 
area;  the  position  of  the  trailing  edge  relative  to  the  wake  center  is 
indicated  on  the  curves.  It  is  seen  from  these  results  that  the  shape 
of  the  turbulent  velocity  profiles  closely  resembles  that  of  the  total- 
head  profiles.  The  main  differences  are  the  more  gradual  decay  of  the 
turbulent  velocities  at  the  outer  edges  of  the  wake  and  the  somewhat 
sharper  peaks  for  the  total-head  loss. 
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Wake  widths  from  figures  3-7  have  been  plotted  as  a  function  of 
angle  of  attack  in  figure  8.  In  defining  wake  width,  the  edge  of  the 

~? 

wake  is  arbitrarily  taken  as  the  point  at  which  total-head  loss  or  w 

is  one-tenth  of  the  maximum  value.  It  is  seen  frcm  figure  8  that  the 

T 
wake  widths  obtained  from  the  total-head  survey  and  the  w  measure- 
ments are  in  close  agreement  for  the  angle-of-attack  range  investigated. 
It  should  be  mentioned  that  wake  widths  based  on  profiles  of  the  RMS 
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turbulent  velocity,  v  w  ,  would  necessarily  be  larger  than  the  present 
wake  widths,  which  are  based  on  the  mean-square  velocity. 

A  comparison  of  the  total-head  profiles  for  angles  of  attack  from 
0  to  lii  degrees  is  presented  in  figure  9.  A  similar  comparison  for 

the  w  measurements  is  given  in  figure  10.   For  both  total  head  and 

2 
w  the  profile  shapes  exhibit  approximate  similarity  for  the  angle- 
of-attack  range  investigated. 

Variation  of  the  wing  drag  coefficient  with  angle  of  attack  is 
shown  in  figure  11.  This  result  was  obtained  by  integration  of  the 
total-head  profiles  measured  one  chord  length  behind  the  wing  trailing 
edge.  The  abrupt  rise  in  drag  in  the  vicinity  of  3  degi^es  angle  of 
attack  apparently  results  from  separation  of  the  boundary  layer  at  the 
leading  edge  of  the  wing.  It  was  observed  during  a  limited  tuft  study 
that  flow  was  extremely  irregular  over  the  entire  upper  surface  of  the 
wing  at  6  degrees  incidence.   Flow  over  the  wing  was  smooth  at  2 
degrees  except  for  a  small  region  of  moderately  irregular  flow  at  the 
leading  edge. 

Variation  of  the  relative  intensity  of  the  downwash  fluctuations 
with  angle  of  attack  is  presented  in  figure  12.  These  values  were 

obtained  by  integrating  the  w  profiles  measured  one  chord  length 
behind  the  wing  trailing  edge.  The  intensity  of  the  downwash  fluc- 
tuations increased  sharply  at  the  higher  angles  of  attack,  probably 
due  to  the  occurrence  of  leading-edge  separation,  as  previously  dis- 
cussed. 

Turbulent  downwash  intensity  is  shown  as  a  function  of  drag 

coefficient  in  figure  13.   Integrated  w  intensity  increased  approxi- 
mately vdth  the  3/2  power  of  the  drag  coefficient  in  these  tests. 
When  the  local  intensity  at  a  given  point  in  the  wake,  such  as  the 
wake  center,  is  plotted  against  drag,  the  variation  is  nearly  linear. 
A  large  part  of  the  increase  in  integrated  intensity  is  due  to  the 
broadening  of  the  wake  at  the  higher  angles  of  attack. 
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Power  Spectra  of  the  Downwash  Fluctuations 


Effect  of  angle  of  attack  and  wake  position.  -  The  resiilts  of 
measurements  of  the  power  spectra  of  the  downwash  fluctuations  one  chord 
length  behind  the  wing  trailing  edge  are  presented  in  figures  li;  to  20 
for  angles  of  attack  from  0  to  li;  degrees.   In  general,  the  power  spectrum 
was  measured  directly  behind  the  wing  trailing  edge  and  also  at  the  outer 
edge  of  the  wing  wake  at  each  angle  of  attack.  The  relative  height  of 
the  turbulent  trace  for  the  two  positions  in  the  wake  may  be  seen  from 
the  wake  profile  of  RMS  downwash  at  the  top  of  each  figure.  The  vertical 
tics  on  the  wake  profiles  show  the  position  at  which  power  spectra  were 
measured. 

A  constant  half -band  width  of  30  cycles  per  second,  based  on  UO 
decibels  attenuation,  was  used  for  all  spectrum  measurements.  The  filter- 
ing characteristics  of  the  wave  analj^er  are  shown  in  figiire  1. 

From  figures  Hx   to  20  it  is  seen  that  wing  angle  of  attack  has  an 
important  effect  on  the  power  spectrum  of  the  downwash  velocity.  The 
spectrum  at  zero  incidence  (fig.  lii.)  is  dominated  by  the  discrete  fre- 
quency corresponding  to  a  vortex  street  shed  from  the  bliuit  trailing 
edge  of  the  wing.  At  3-1/2  degrees  angle  of  attack,  the  trailing  edge 
vortex  street  is  still  apparent,  but  the  bulk  of  the  txirbulent  energy  is 
spread  over  a  broad  band  of  frequencies. 

With  further  increase  in  angle  of  attack,  the  downwash  spectrum 
near  the  wake  center  shows  a  progressive  decrease  in  the  predominant  fre- 
quencies. The  bulk  of  the  turbulent  energy  is  concentrated  near  the 
wake  center,  and  so  a  general  conclusion  from  this  data  is  that  wake  fre- 
quency decreases  as  wing  incidence  is  increased. 

By  lij.  degrees  angle  of  attack  (fig.  20),  a  vortex  street  is  formed 
behind  the  stalled  wing,  and  the  measured  spectrum  corresponds  essen- 
tially to  a  delta  function  which  is  spread  out  by  the  finite  band  width 
of  the  wave  analyzer  (fig.  1).  At  first  thought,  it  seems  imusual  that 
the  frequency  measured  near  the  center  of  the  vortex  street  is  not  double 
the  frequency  at  the  wake  edges  because  of  the  influence  of  both  rows  of 
the  vortex  street.  Therefore,  a  calculation  of  the  velocities  induced 
by  a  vortex  street  was  made  using  the  equations  of  reference  lit.  Thei 
results  are  presented  in  figure  21,  along  with  the  geometry  of  the  as- 
sumed vortex  street.  The  upper  graph  in  figure  21  gives  the  variation 
of  the  longitudinal  induced  velocity,  u,  with  distance  along  the  vortex 
street,  and  the  lower  graph  presents  the  variation  of  induced  downwash, 
-w.  The  induced  velocities  were  calculated  along  the  center  of  the  vor- 
tex street  (z  =  0)  and  also  along  a  line  located  above  the  vortex  street 
(z  =  b).  As  expected,  the  predominant  frequency  of  the  u  velocity  in 
the  center  of  the  wake  is  double  the  frequency  outside  the  vortex  street. 
However,  the  downwash  velocity  in  the  wake  center  does  not  show  fre- 
quency doubling;  this  result  is  in  agreement  with  the  measurements  of 
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the  downwash  power  spectrum  presented  in  figure  20. 

Below  the  stalling  incidence,  the  downwash  spectrum  is  much 
broader  near  the  center  of  the  wake  than  at  the  outer  edge  of  the 
wake.  Moreover,  the  predominant  frequencies  at  the  edge  of  the  wake 
are  practically  independent  of  wing  incidence  for  angles  of  attack 
from  6  to  lU  degrees.  Because  there  is  very  little  turbulent  energy 
in  the  outer  edges  of  the  wake,  subsequent  analysis  of  the  data  will 
be,  for  the  most  part,  based  on  the  spectrum  measurements  directly 
behind  the  wing  trailing  edge  (2=0). 

Relation  between  frequency  and  wake  width.  The  data  presented 
in  figures  lU  to  20  show  that  the  predominant  frequencies  in  the  wing 
wake  generally  decrease  as  the  wake  width  increases.  In  order  to 
illustrate  this  result  more  explicitly,  it  is  convenient  to  define 
an  average  frequency  for  the  turbulent  power  spectrum.  The  average 
wake  frequency,  n  ,  is  defined  as  the  frequency  which  divides  a 

power  spectnmi  into  two  equal  areas.   In  figure  22  the  reciprocal  of 
the  average  wake  frequency  is  plotted  against  an  effective  wake  width, 
defined  as  the  wake  width  divided  by  average  wake  velocity.  The  ex- 
periments show  that  the  average  wake  frequency  is  inversely  propor- 
tional to  the  effective  wake  width.  In  other  words,  a  non-dimensional 
frequency  or  Strouhal  number  based  on  wake  width  and  wake  velocity  is 
nearly  constant.  This  result  is  in  agreement  with  a  recent  correla- 
tion based  on  free-streamline  theory  for  a  series  of  bluff  bodies 
(ref.  1$). 

A  line  labeled  "instability  theory"  in  figure  22  was  calculated 
visjjig  the  instability  theory  of  wakes  presented  in  reference  I6.  The 
possibility  of  applying  stability  theory  to  a  turbulent  velocity 

2 

profile  has  never  been  fiilly  discussed.   Instability  of  laminar  shear 

flows  generally  implies  the  eventual  formation  of  a  turbulent  flow. 
When  the  flow  is  already  turbulent,  the  meaning  of  an  instability 
calculation  is  far  from  clear.  However,  even  in  laminar  flow  only  the 
very  beginning  of  instability  can  be  handled  by  a  linear  theory;  it 
is  not  impossible  that  the  application  of  similar  considerations  to 
turbulent  flow  could  provide  interesting  results.  For  example,  it  is 
known  that  Kanr.an  vortex  streets  and  three-dimensional  Taylor  vortices 
may  exist  even  in  turbulent  flow.  Therefore,  it  is  interesting  to 
compute  the  frequencies  predicted  for  neutral  stability  and  to  com- 
pare these  frequencies  with  the  measured  results,  as  in  figure  22. 
Instability  theory  predicts  the  same  inverse  relationship  between 
frequency  and  effective  wake  width  shown  by  the  experimental  measure- 
ments. Although  the  actual  magnitude  of  the  predicted  frequency  is 


^ = 

A  few  remarks  concerning  the  existence  of  a  large-scale   structure 
iji  turbulent  shear  flows  have  been  made  in  reference  13 . 


10  NACA  TM  11+27 


in  close  agreement  with  the  measured  results,  this  is  very  likely  coin- 
cidental. First  of  all,  the  theory  assumes  neutrally  stable  disturb- 
ances whereas  amplified  disturbances  woxild  probably  be  more  important 
in  any  real  flow.  Also,  different  average  frequencies  coiild  be  ob- 
tained from  the  measiored  spectra  depending  on  the  method  chosen  to 
define  an  average.  Details  of  the  stability  calculation  are  presented 
in  Appendix  B. 

rtelation  between  frequency  and  projected  wing  chord.  -  It  has  been 
observed  in  several  investigations  that  the  Strouhal  number  based  on 
the  projected  chord  of  a  wing  is  nearly  constant  for  angles  of  attack 
above  the  stall  (ref.  1).  In  order  to  deteimine  whether  such  a  constant 
Strouhal  number  is  obtained  for  smaller  angles  of  attack,  the  average 
wake  frequencies,  n  ,  obtained  from  measured  power  spectra  are  shown 

in  figure  23  as  a  function  of  projected  airfoil  breadth.  This  projected 

breadth  is  defined  as  the  total  airfoil  projection,  which  includes  the 

finite  thickness  of  the  wing.  The  line  labeled  "empirical"  corresponds 

to  a  constant  projected  Strouhal  number  obtained  by  extrapolating  the 

frequency  measured  for  the  stalled  wing  (  a  =  lli  ). 

The  empirical  extrapolation  gives  a  reasonable  prediction  of  the 
average  wake  frequency  except  for  an  angle  of  attack  of  3-1/2  degrees. 
At  this  angle  of  attack  the  airfoil  does  not  have  the  leading-edge 
separation  that  occurred  at  higher  angles,  and  so  it  is  not  surprising 
that  wake  frequency  is  not  related  to  projected  chord.  When  frequency 
is  related  to  wake  width  instead  of  projected  chord,  the  unseparated 
flow  condition  does  not  present  this  difficulty.  Therefore,  it  appears 
that  wake  width  is  more  suitable  than  projected  chord  for  correlating 
frequency  measurements. 


Power  Spectra  of  the  Longitudinal  and  Spanwise  Velocities 


Although  downwash  fluctuations  are  generally  more  effective  as  a 
cause  of  buffeting  loads,  the  power  spectra  of  the  longitudinal  and 
spanwise  velocities  are  necessary  for  more  complete  understanding  of 
turbulent  wakes.  Longitudinal  velocity  spectra  were  measured  for  an 
angle  of  attack  of  6  degrees  at  several  vertical  positions  in  the  wake, 
and  the  results  are  shown  in  figures  2k   and  25.  The  corresponding 
downwash  spectra  are  presented  in  figiures  26  and  27. 

2      2 

Looking  at  the  data  as  a  whole,  the  u  and  w  power  spectra  gen- 

nc 
erally  have  pronoimced  maxima  at  values  of  tj^  =  0,65  or  1.30.  The 

2      2  °° 

u  and  w  spectra  are  similar  to  each  other  over  the  outer  portion  of 

the  wake,  but  the  spectra  are  markedly  different  near  the  wake  center. 
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In  order  to  determine  whether  the  wing  wake  resembles  a  vortex 
street  at  this  angle  of  attack,  the  velocities  induced  by  a  vortex 
street  (fig.  21)  were  expressed  in  terms  of  their  Foiirier  coefficients. 
The  resulting  hairaonic  analysis  of  an  idealized  vortex  street  is 
presented  in  figure  28  for  several  vertical  locations  above  the  center 
of  the  vortex  street.  The  first  and  second  Fourier  coefficients  for 
the  u  velocity  are  on  the  left  of  the  figure,  and  the  corresponding 
downwash  coefficients  are  on  the  right.  Comparing  the  calculated 
results  with  the  measurements  given  in  figures  2k   to  27,  it  is  seen 
that  the  measured  data  do  not  show  consistent  agreement  with  the  cal- 
culations. Therefore,  it  seems  likely  that  the  structure  of  the  wing 
wake  at  low  angles  of  attack  is  more  complicated  than  that  of  a  simple 
vortex  street. 

2 

In  figure  29  the  u  spectrum  directly  behind  the  wing  trailing 

edge  is  compared  with  the  spectrum  of  the  spanwise  velocity  component. 
These  spectra  may  be  compared  with  the  spectra  for  isotropic  turbu- 

2 

lence.  The  longitudinal  or  u  spectrim  of  isotropic  turbulence  has 

been  obtained  experimentally  in  the  investigations  of  references  6 
and  17.   From  such  resiilts  the  transverse  spectra  may  be  calculated 
using  the  continuity  relation  for  isotropic  turbulence  (ref.  18). 
Measurements  of  the  longitudinal  spectrum  may  be  approximated  by  the 
relation 


F-,(n)  =  ^  — ^  (2) 

^       ^00  (1  +  C) 


and  the  corresponding  transverse  spectrum  is 


00  {1  +  K  ) 


The  variable  C  is  a  non-dimensional  frequency  given  by 

00 

and  L  is  the  turbulent  macros cale  of  the  longitudinal  velocity 

U  F^ (0) 
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The  isotropic  power  spectra  shown  in  figure  30  were  calculated 
for  the  same  value  of  L  as  that  measured  in  the  wing  wake,  figure  29. 

2      2 

Comparison  of  figures  29  and  30  indicates  that  the  u  and  v  spectra 

near  the  center  of  the  wing  wake  are  similar,  respectively,  to  the 
longitudinal  and  transverse  spectra  of  isotropic  turbulence.  This 
similarity  was  not  shown  at  other  vertical  positions  in  the  wake,  al- 
though part  of  the  total  turbulent  energy  may  still  be  attributable  to 
an  isotropic-like  turbulence. 


Effect  of  Free-Stream  Velocity 


The  results  presented  so  far  were  obtained  at  a  free-stream 
velocity  of  about  50  feet  per  second.  In  order  to  determine  whether 
the  results  may  be  used  for  other  free-stream  velocities,  representa- 
tive tests  were  made  at  a  lower  velocity.  These  additional  tests  do 
not  provide  significant  information  on  the  effect  of  Reynolds  nianber 
variation,  but  should  facilitate  interpretation  of  the  results.  For 
example,  it  is  generally  assumed  that  aerodynamic  frequencies  are 
proportional  to  velocity,  and  so  this  assimption  should  be  checked. 

Wake-energy  distribution.  -  The  fluctuating-downwash  profile  is 
shown  in  figure  31  for  velocities  of  30  and  k9   feet  per  second.  Although 
the  profile  is  somewhat  steeper  for  the  lower  velocity,  the  wake  width 
is  practically  unaffected  by  the  velocity  change.  The  overall  effect  of 
the  velocity  change  on  the  profile  is  negligible. 

Discrete  frequencies.  -  The  effect  of  velocity  on  the  periodic 
walce  beliind  the  blunt  trailing  edge  of  the  flat-plate  airfoil  is  shovm 
in   figure  32  for  zero  angle  of  attack.  At  the  higher  velocities,  fre- 

quency  is  proportional  to  velocity.  The  Strouhal  number,  -rr—  ,  of  the 

CO 

fluctuations,  based  on  trailing-edge  thickness,  is  0.2l|,  which  is  in 
good  agreement  viith  the  Strouhal  number  of  0.23  obtained  in  reference  19 
from  schlieren  observations  of  the  vortex  street  behind  a  blunt  trailing- 
edge  airfoil.  For  velocities  below  about  20  feet  per  second  (fig.  32), 
the  frequencies  in  the  wing  wake  become  nonlinear  with  velocity.  This 
corresponds  to  a  decrease  in  Strouhal  nimber  for  Reynolds  nimbers 

below  10  ,  with  Reynolds  number  based  on  trailing-edge  thickness.  The 
Strouhal  number  of  .a  cylinder  also  shows  a  decrease  for  Reynolds  num- 
bers below  10^.   (See  fig.  li;9,  ref.  20.) 

Velocity  fluctuations  in  the  wake  were  also  periodic  above  the 
stalling  angle  of  the  wing.  The  effect  of  free-stream  velocity  on  wake 
frequency  is  shown  in  figure  33  for  two  angles  of  attack.  Over  the 
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Reynolds  nimber  range  of  these  tests,  frequency  is  proportional  to 
velocity.  The  Strouhal  nximber  based  on  projected  chord  is  0.19  for 

a  =  1U.5°  and  O.I8  for  a  =  21.7°.   In  the  discussion  of  the  accuracy 
of  measurements,  it  was  noted  that  these  frequencies  are  probably 
slightly  high  due  to  wall  interference.  Nevertheless,  the  values  are 
well  within  the  range  of  Strouhal  numbers  obtained  in  various  inves- 
tigations of  vortex  streets  behind  stalled  plates  and  airfoils  (ref.  1), 

Power  spectra.  -  The  power  spectrum  of  downwash  fluctuations 

behind  the  wing  at  a  moderate  angle  of  attack  is  shown  in  figure  3U 

for  two  values  of  the  free-stream  velocity.  Roughly  speaking,  the 

normalized  power-spectral  density  is  uniquely  dependent  on  the  usual 

nc 
frequency  parameter,  tj—  ,  for  the  two  velocities.  Some  of  the  dis- 

ao 

agreement  between  the  two  spectra  very  likely  results  from  the  replace- 
ment of  hot-wires  during  the  interval  between  the  two  tests.  In  any 
case,  the  overall  agreement  of  the  spectra  appears  quite  satisfactory. 


DISCUSSION 
Prediction  of  RMS  Lift  in  Tail  Buffeting 


Formulas  for  the  prediction  of  the  RMS  lift  fluctuations  on  a 
tail  surface  located  in  a  turbulent  wake  have  been  given  in  reference 
21.  However,  in  order  to  obtain  explicit  results,  it  was  necessary 
to  assume  that  the  turbulence  striking  the  tail  surface  was  isotropic. 
The  effect  of  this  assumption  may  be  evaluated  by  using  an  actual 
power  spectrum  measured  in  these  tests  to  calculate  RMS  lift  on  a 
hypothetical  tail  surface. 

A  convenient  relation  for  calculating  mean-square  lift  was  de- 
rived in  reference  21,  and,  except  for  a  slight  difference  in  notation, 
the  relation  is: 


00 
C^^  =  h   nV  JFCn)  y   (k)|^  dn 


(6) 


The  assumptions  on  which  this  equation  is  based  are:   (1)  the  tail 
aerodynamics  may  be  described  by  a  two-dimensional  admittance  func- 
tion, ^{k);    (2)  the  scale  of  turbulence  is  large  compared  with  the 
tail  span. 

It  is  evident  from  equation  6  that  certajji  parameters,  such  as 


2h 
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velocity  and  wing  chord,  must  be  assumed  in  order  to  calculate  the 
mean-square  lift  coefficient.  However,  the  ndnimum  number  of  assump- 
tions are  made  if  equation  6  is  modified  so  as  to  incorporate  a  useful 
similarity  property  of  the  power  spectrumj  namely,  that 


U  F(n) 

00 


(7) 


as  seen  from  figure  3U.  Then  equation  6  becomes 


1  2  2 


(8) 


CO' 


The  ob^iect  of  the  lift  calculation  may  now  be  described  more 
explicitly.  Using  the  same  tail  surface  and  aerodynamic  admittance 


function, 


Cj^  /  a 


is  to  be  calculated  for  a  power  spectrum  measured 


behind  the  wing  (fig.  17,  z  =  O)  and  compared  vjith  a  similar  result 
calculated  for  the  transverse  isotropic  spectrum  (fig.  30).  The 
aerodynamic  admittance  of  a  tail  having  c./c  =0.2  was  calculated  using 

equation  13  of  reference  3,  and  the  result  is  shown  in  figure  35-  For 
this  tail  surface,  the  comparison  of  RMS  lift  for  an  actual  wing  down- 
wash  spectrum  vdth  the  B11S   lift  for  isotropic  turbulence  is  as  follows: 


SPECTRUM 

VJing 
Isotropic 

0.37 
O.Ul 
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I^Vom  this  comparison,  it  is  evident  that  the  exact  shape  of  the 
turbulent  power  spectrum  has  a  relatively  small  effect  on  the  mean 
lift  fluctuations  experienced  by  a  tail  siirface.  A  simple  closed-form 

2  ~2 

expression  is  available  for  calculating  C,  /a    in  isotropic  tur- 

bulence  (equation  28  of  reference  21),  and  the  above  comparison  shows 

that  this  equation  should  be  quite  xiseful  for  predicting  the  RMS  lift 

on  a  tail  in  the  wake  of  a  wing.  However,  in  order  to  apply  the 

formula,  it  is  necessary  to  estimate  the  macroscale  of  the  turbulence 

striking  the  tail  surface.  For  the  comparison  just  presented,  the 

2 
turbulent  scale  vas  obtained  from  a  measured  u  spectrum  using  equation 

5.  Such  results  are  not  usually  available  to  designers,  and  so  it  is 

necessary  to  have  a  simpler  means  of  estimating  the  scale  of  tiirbulence 

in  a  wing  wake.  A  possible  method  for  obtaining  such  an  estimate  will 

now  be  discussed. 

The  turbulent  scale,  L,  may  be  used  to  define  a  characteristic 
frequency  by  simple  dimensional  analysis: 

U 
n=  -22.  (9) 

2nL 

It  might  be  expected  that  this  frequency  would  correspond  to  the  pre- 
dominant frequency  in  the  turbulence,  and  this  assumption  may  be 
checked  using  the  data  of  the  present  tests.  The  characteristic  fre- 

2 

quency  calculated  from  the  u  spectrum  of  figvire  29  and  equation  9 

was  found  to  be  nc/U   =  1.2U.  This  frequency  corresponds  almost 

00  2 

exactly  to  the  predominant  frequency  in  the  w  power  spectrum,  figure 
17.  Therefore,  equation  9  provides  a  simple  means  of  estimating  tur- 
bulent scale  because  wake  frequencies  can  be  obtained  more  easily 
than  a  power  spectrvim.  For  example,  a  tuft  probe  could  be  placed  in 
the  wing  wake  and  the  mean  frequency  could  then  be  obtained  from 
motion  pictures  of  the  tuft  movement. 


Buffeting  Measurements  with  Detector  Airfoils 


From  the  results  of  some  earlier  investigations  of  tail  buf- 
feting (ref.  22,  23),  it  was  concluded  that  lift  fluctuations  may  be 
obtained  vrell  outside  the  limits  of  the  total-head  wake.  This  res\ilt 
appears  to  be  at  variance  vri.th  the  present  investigation  and  that  of 
reference  21;.  Therefore,  it  seems  worthi-jhile  to  discuss  the  reasons 
for  the  different  conclusions. 

The  fact  that  relatively  strong  lift  fluctuations  have  been 
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detected  beyond  the  region  of  the  total-head  wake  is  not  surprising 
for  several  reasons.  First,  the  intensity  of  the  fluctuations  observed 
in  references  22  and  23  was  taken,  roughly  speaking,  as  the  maximuni 
amplitude  of  the  detector  airfoils.  Therefore,  the  data  obtained  in 
this  manner  represent  an  upper  limit  of  the  lift  fluctuations  and  do 
not  necessarily  represent  a  statistical  average  of  the  fluctuations. 
In  the  present  investigation,  it  was  generally  observed  that  near  the 
edge  of  the  total-head  wake,  occasional  bursts  of  turbulence  were  found 
to  be  notably  larger  than  the  general  level  of  the  turbulent  velocity 
fluctuations.  l-Jhen  such  bursts  occur  in  the  outer  edges  of  the  wake, 
a  wake  vri.dth  based  on  maximum  values  would  be  larger  than  the  width 
obtained  from  statistically  averaged  data. 

Even  in  the  absence  of  such  bursts,  the  profiles  in  references 
22  and  23  would  be  broader  than  those  of  the  present  report  as  only 
the  latter  are  expressed  in  terras  of  mean-square  fluctuations.  The 
profiles  of  references  22  and  23  correspond  more  nearly  to  RMS  fluctua- 
tions, and  RMS  values  must  necessarily  fall  off  more  rapidly  than  mean 
squares. 

Finally  it  is  noted  that  the  predominant  velocity  fluctuations 
near  the  center  of  the  wake  occur  at  higher  frequencies  than  those 
observed  at  the  t-rake  edges  (figs.  26  and  27).  Because  low  frequency 
velocity  fluctuations  are  more  effectively  converted  to  lift  fluctua- 
tions (ref.  3),  the  lower  frequencies  found  in  the  outer  portions  of 
a  wake  are  less  attenuated  than  the  high  frequency  energy  in  the  wake 
center.  This  effect  supplies  another  contributing  cause  for  the  breadth 
of  the  fluctuating  lift  profiles  presented  in  references  22  and  23. 

In  order  to  illustrate  the  effects  just  discussed,  representative 
data  obtained  in  the  present  investigation  have  been  shown  in  several 
different  forms  in  figures  36  to  38.  Using  the  measured  power  spectra 
for  the  dovjnwash  fluctuations,  figure  36,  the  frequency  distribution 
of  RMS  lift  has  been  estimated  for  an  airfoil  located  at  various 
heights  in  the  wing  wake.  The  simulated  horizontal  tail  had  a  chord 
equal  to  20  percent  of  the  wing  chord.  The  form  of  the  approximate 
admittance  function  used  to  calculate  lift  spectra  is  shown  in  figure 
35.  On  the  basis  of  results  from  reference  k,   it  is  expected  that 
such  an  admittance  function  is  representative  of  the  admittance  of 
the  sensitive  airfoils  used  by  Ferri  (ref.  23).  The  resulting  fre- 
quency distribution  of  RIIS  lift  estimated  in  this  manner  is  presented 
in  figure  38.  A  comparison  of  the  relative  intensities  obtained  from 
the  RMS  lift  spectra  (fig.  38)  and  the  conventional  angularity  spectra 
(fig.  36)  is  as  folloTS: 
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From  this  result,  it  is  seen  that  tests  using  sensitive  airfoils  would 
be  expected  to  indicate  greater  wake  widths  than  those  obtained  from 
hot-vire  siirveys.  The  occurrence  of  bursts,  as  previously  disciissed, 
could  cause  even  greater  differences  betv;een  lift  and  angularity  data. 


CONCLUSIONS 


1.  The  mean  square  of  the  dovmwash  fluctuations  v^as  found  to 
have  approximately  the  same  variation  with  vertical  position  in  the 
wing  wake  as  that  shown  by  measurements  of  total -head  loss. 

2.  The  intensity  of  downwash  fluctuations  increased  rapidly  when 
the  flow  separated  from  the  wing  surface.  For  the  present  tests, 
integrated  do-tmwash  intensity  increased  with  approximately  the  3/2 
power  of  idng  drag  coefficient. 

3.  Power  spectra  for  the  downvrash  velocity  were  generally  broader 
at  the  low  angles  of  attack  than  at  high  angles.  The  average  fre- 
quency in  the  wake  was  proportional  to  the  ratio  of  wake  velocity  to 
T-rake  ladth.  Vake  vridth  appears  to  be  a  more  suitable  parameter  for 
correlating  spectrum  measurements  than  the  projected  airfoil  chord, 

U.  The  outer  edges  of  the  wing  wake  contained  relatively  little 
high-frequency  energy,  idiereas  the  turbulent  energy  near  the  wake 
center  was  distributed  over  a  fairly  vdde  range  of  frequencies. 

5.  The  wake  structure  at  low  angles  of  attack  appears  to  be 
much  more  complicated  than  for  the  stalled  airfoil,  for  which  a  vor- 
tex street  forms  behind  the  viing, 

6.  Po>7er  spectra  measured  at  different  velocities  were  found 

to  be  similar  vhen   frequency  and  povjer  spectral  density  were  non- 

dimensionalized  by  the  ratio  of  airfoil  chord  to  free-stream  velocity, 

c/U  . 
00 
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7.  On  the  basis  of  calculated  results,  the  exact  shape  of  the 
turbulent  poi^rer  spectrum  was  found  to  have  a  relatively  small  effect 
on  the  RMS  lift  of  a  tail  surface  located  in  the  vra.ke  of  the  wing. 
Formulas  for  predicting  lift  assuming  isotropic  turbulence  should  give 
satisfactory  results  if  the  turbulent  scale  is  correctly  estimated. 


APPENDIX  A 

Syitibols 

The  following  symbols  are  defined  as  they  are  used  in  the  main 
body  of  this  paper.  Certain  letters  have  been  used  with  different 
meaning  vrithin  Appendix  B  in  order  to  retain  the  notation  of  reference 
16. 

a  distance  between  vortices  in  the  same  row 

b         breadth  of  vortex  street 

C^        drag  coefficient,  ^^^gAP^ 
D  ^  *    qc 

C,         lift  coefficient,  ^^^V^P^" 
L  '    qc 

c         airfoil  chord 

F(n)       power  spectral  density  of  a  random  variable j 


00 

F(n)  =  U  /  R(t)cos  2nnT  dx 
'o 


'"'  ■  'I 


F^(n)      power  spectral  density  for  longitudinal  velocity  in  isotro- 
pic turbulence 

FpCn)      power  spectral  density  for  transverse  velocity  in  isotro- 
pic turbulence 

k         reduced  frequency,  coc/2U 

L         turbiilent  macroscale,  equation  5 
n         frequency,  cycles/second 
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n-       the  first  harmonic  frequency  of  a  vortex  street 

n        average  vake   frequency,  defined  as  the  frequency  vxhich 
divides  a  power  spectrum  into  two  equal  areas 

q        dynamic  pressure,  -k  pU 

R(t:)      time-correlation  function  of  a  random  variable; 

R(t)  =  u(t)  u(t  +  t)/  u^ 

pOO 

R(t^)  =  /    F(n)  cos  2nnT  dn 


TE       trailing  edge 
t        m.ng   thickness 


U        free-stream  velocity 
00 


u  longitudinal  perturbation  velocity 

V  spanwise  perturbation  velocity 

w  vertical  perturbation  velocity 

W  wake  i-jidth 

X  distance  along  stream  direction  from  vring  trailing  edge 

z  vertical  distance  above  centerline  of  vdng  trailing  edge 

a  angle  of  attack 

/""  vortex  strength 

Ah  difference  betvreen  free-stream  and  local  total  head 

Ly  height  above  wake  center 

K  non-dimensional  frequency,  oiL/U 

p  free  stream  density 
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T        time 

^(k)     ratio  of  lift  at  frequency  k  to  lift  at  zero  frequency 

CO        circular  frequency,  2nn 

APPENDIX  B 


Application  of  Stability  Theory  to  the  Prediction  of 
Wake  Frequencies 


In  this  appendix  assumptions  underlying  the  stability  theory  of 
vakes  (ref.  16)  are  briefly  outlined.  Also,  certain  additional  formulas 
are  developed  in  order  to  facilitate  comparison  between  theory  and 
experiment. 

Hollingdale  (ref.  16)  considers  a  mean  wake  flow  ■vAose  velocity 
relative  to  the  stream  at  infinity  is 


(Bl) 


With  the  assumption  of  small  two-dimensional  disturbances  a  linear, 
homogeneous  equation  for  the  disturbance  velocities  is  obtained  from 
the  incompressible  Navier-Stokes  equations.  The  stream  function  \|f 
may  be  assumed  periodic,  so  that 

t  =  ,,((z)e^"^^-^*^  (B2) 

The  wave  length  of  the  disturbance  is  2Ti/a,  and  c  is  the  wave  velocity. 
The  basic  equation  for  determining  the  stability  of  small  distiorbances 
in  a  wake  or  other  region  of  shear  flow  is  then 

(U  -  c){^»   -  a^^)  -   UV  =  ^  (j^f'^  _  2  aV"  +  a^)   (B3) 

From  the  complete  stability  equation  B3  Hollingdale  has  deter- 
mined that  the  wave  velocity  c  for  an  undamped  distiu-bance  is  equal 
to  the  value  of  U  at  the  point  of  inflection  of  the  velocity  profile 
U(z).  This  resiilt  is  valid  for  symmetrical  or  monotonic  profiles  in 
which  U  approaches  a  constant  for  large  z. 
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For  sufficiently  large  Reynolds  nuinbers  the  viscous  terms  of 
equation  B3  rnay  be  omitted: 

(U  -  c)((^»  -  aV)  -  UV  =  0  (BU) 


The  velocity  profile  may  be  assumed  to  be  of  the  form 


U  =  A(l  +  cos  kz)  0  $  I  z  I  ^   nA 

U  =  0  |z| >  n/k 


(B5) 


v^ere  n/k  is  the  wake  semiwidth,  subsequently  denoted  as  W/2.  The 
constant  A  must  equal  the  wave  velocity  in  order  that  U  =  c  at  the 
point  of  inflection  U"  =  0.   For  the  particiilar  form  of  profile 
chosen  the  stream  function  must  therefore  satisfy  the  relations 

j!^"  +  (k^  -  a^)jz{  =  0        0$  |z  I  ^   nA 

P  (B6) 

^»   _  aV  =  0  |z|  >  nA 


The  boundary  conditions  to  be  satisfied  by  ^   are  that  it  is  an  even 

function  in  the  region  (  -  rz   j  tt  )  and  that  it  approaches  zero  for 

large  z.  Solutions  satisfying  these  conditions  and  the  differential 
equations  B6  are 


jzf  =  P  cos  yjk^  -  a^     z  0  $  I  z  I    ^      nA 

^  =  Qe"''l^l  |z|>  nA 


These  equations  provide  the  characteristic  equation  for  determining 
the  wave-length  parameter  a  if  /  and  ^'  are  required  to  be  continuous 
at  Izl  =  nAi 


tan    riVk     -a  ^    _^^^  ^^g^ 
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The  parameter 

_     _     a     _       «  (3^) 

k 


is  simply  the  ratio  of  the  wake  tiddth  to  the  wave  length  of  the  disturb- 
ance.  The  important  result  from  eqxiation  B8  is  that  this  ratio  is  a 
constant  and,  in  fact,  has  the  numerical  value 


a     =  0.926  (BIO) 


Disturbance  waves  of  the  type  predicted  by  Hollingdale  reach  a 
point  v;hich  is  fixed  relative  to  the  body  producing  the  wake  with  the 
frequency 


n  =  ^  (U^  -  A)  (BID 

dn         00 


where  U   is  the  free-stream  velocity  and  A  is  the  wave  velocity.  In- 
troducing Hollingdale ' s  result  for  the  proportionality  between  wave 
length  and  wake  width,  equation  B9,  provides  the  expression  for  the 
frequency  of  small  periodic  disturbances: 


(U^  -  A)  a 
IT 


n=-22_ : (B12) 


Although  the  value  of  the  wave  velocity  A  may  be  determined  as 
the  velocity  at  the  point  of  inflection  of  the  velocity  profile,  it  is 
difficult  in  practice  to  determine  the  inflection  point  of  an  experi- 
mental velocity  profile.  Hence,  it  is  desirable  to  have  a  more  practi- 
cal method  of  estimating  wave  velocity.  The  drag  may  be  obtained  by 
integration  of  the  total-head  loss  if  the  survey  plane  is  a  sufficient 
distance  from  the  wing  trailing  edge: 


C^  =  .  AH 


(I)  («") 
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The  relation  between  total-head  loss  and  the  velocity  deficit  U  in  the 
wake  may  be  obtained  from  Bernoulli's  equation  by  assuming  the  static 
pressure  in  the  wake  is  equal  to  the  free-stream  pressure. 


Ah  _  u   /„   u  \  (g^) 


When  equation  B5  for  the  velocity  variation,  is  substituted  into 
equations  B13  and  Blii,  a  relation  for  the  wave  velocity  in  terms  of 
drag  and  wake  width  is  obtained: 


■-F 


3  A    _  ,    L       3    S° 


?  U      ^  -  V""  "?  IT 

00 


(B15) 
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Figure  1.-  Filtering  characteristics  of  Hewlett-Packard  model  5OOA  wave 
analyzer  for  50  cps  half -band  width. 
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Figure  lk.~   RMS  downwash  intensity  and  w   power  spectrum  measured 
1  chord  length  behind  trailing  edge;  a  =  0;  V^  =  ^2   fps. 
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Figure  15 •-  RMS  downwash  intensity  and  m       power  spectrum  measured 
1  chord  length  behind  trailing  edge;  a  =  3°30 ' ,   Uoo  =  52  fps . 
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Figure  l6.-  RMS  downwash  intensity  and  w'' 
1  chord  length  behind  trailing  edge;  a 


power  spectrum  measured 
5015',  Uoo  =  52  fps. 
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Figure  17- -  RMS  downwash  intensity  and  w^ 
1  chord  length  behind  trailing  edge;  a 


power  spectrum  measured 


6°20' 


Uoo  =  5^  fps. 
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Figure  l8.-  RMS  downwash  intensity  and  w   power  spectrum  measured 
1  chord  length  behind  trailing  edge;  a  =  8°35 ' ;   U^^  =  52  fps. 
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Figure  19--  RMS  downwash  intensity  and  w   power  spectrum  measured 
1  chord  length  behind  trailing  edge;  a  =  10°15 ' ,   Uoo  =  57  fps  for 
z  =  0;  a  =  10°35S   Uoo  =  ^^-8  fps  for  z/c  =  -O.h. 
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Figiire  20.-  RMS  downwash  intensity  and     w^ 
1  chord  length  behind  trailing  edgej   a 
Uoo  =  52  fps  for  z/c   =  -O.k. 
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Figure  21.-  The  velocities  induced  by  a  vortex  street. 
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Figure  22.-  Dependence  of  average  wake  frequency  on  effective  wake  width. 
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Figure   2k.-  v.       power   spectrum;   a  =  6°20 ' ,      Uoo  =  5^  fps,      x   =  Ic 
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Figure  25--  u       power   spectrixm;   a  =  6°20 ' ,      Uoo  =  5^  fps,      x  =  Ic. 
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Figure  26.-  w       power  spectrum;   a  =  6°20 ' ,      Uoo  =  5^  fps. 
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Figure   27.-  w       power  spectrum;   a  =  6°20 ' ,      Uoo  =  5^1  fps,      x  =  Ic. 
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